Superkicks in ultrarelativistic encounters of spinning black holes

We study ultrarelativistic encounters of two spinning, equal-mass black holes through simulations in full numerical relativity. Two initial data sequences are studied in detail: one that leads to scattering and one that leads to a grazing collision and merger. In all cases, the initial black hole spins lie in the orbital plane, a configuration that leads to the so-called "superkicks". In astrophysical, quasicircular inspirals, such kicks can be as large as ~3,000 km/s; here, we find configurations that exceed ~15,000 km/s. We find that the maximum recoil is to a good approximation proportional to the total amount of energy radiated in gravitational waves, but largely independent of whether a merger occurs or not. This shows that the mechanism predominantly responsible for the superkick is not related to merger dynamics. Rather, a consistent explanation is that the "bobbing" motion of the orbit causes an asymmetric beaming of the radiation produced by the in-plane orbital motion of the binary, and the net asymmetry is balanced by a recoil. We use our results to formulate some conjectures on the ultimate kick achievable in any black hole encounter.Comment: 10 pages, 6 figures, 2 table

The BSSN formulation is a partially constrained evolution system

Relativistic simulations in 3+1 dimensions typically monitor the Hamiltonian and momentum constraints during evolution, with significant violations of these constraints indicating the presence of instabilities. In this paper we rewrite the momentum constraints as first-order evolution equations, and show that the popular BSSN formulation of the Einstein equations explicitly uses the momentum constraints as evolution equations. We conjecture that this feature is a key reason for the relative success of the BSSN formulation in numerical relativity.Comment: 8 pages, minor grammatical correction

Can Black Holes Decay to Naked Singularities?

We investigate thermodynamic properties of two types of asymptotically anti-de Sitter spacetimes: black holes and singular scalar field spacetimes. We describe the possibility that thermodynamic phase transitions can transform one spacetime into another, suggesting that black holes can radiate to naked singularities.Comment: 5 pages, Essay for 2001 Gravity Research Foundation competition, to appear in IJMP

Interpolating compact binary waveforms using the singular value decomposition

Compact binary systems with total masses between tens and hundreds of solar masses will produce gravitational waves during their merger phase that are detectable by second-generation ground-based gravitational-wave detectors. In order to model the gravitational waveform of the merger epoch of compact binary coalescence, the full Einstein equations must be solved numerically for the entire mass and spin parameter space. However, this is computationally expensive. Several models have been proposed to interpolate the results of numerical relativity simulations. In this paper we propose a numerical interpolation scheme that stems from the singular value decomposition. This algorithm shows promise in allowing one to construct arbitrary waveforms within a certain parameter space given a sufficient density of numerical simulations covering the same parameter space. We also investigate how similar approaches could be used to interpolate waveforms in the context of parameter estimation.Comment: 5 pages, 3 figures, presented at the joint 9th Edoardo Amaldi Conference on Gravitational Waves and 2011 Numerical Relativity - Data Analysis (NRDA) meetin

Comparisons of binary black hole merger waveforms

This a particularly exciting time for gravitational wave physics. Ground-based gravitational wave detectors are now operating at a sensitivity such that gravitational radiation may soon be directly detected, and recently several groups have independently made significant breakthroughs that have finally enabled numerical relativists to solve the Einstein field equations for coalescing black-hole binaries, a key source of gravitational radiation. The numerical relativity community is now in the position to begin providing simulated merger waveforms for use by the data analysis community, and it is therefore very important that we provide ways to validate the results produced by various numerical approaches. Here, we present a simple comparison of the waveforms produced by two very different, but equally successful approaches--the generalized harmonic gauge and the moving puncture methods. We compare waveforms of equal-mass black hole mergers with minimal or vanishing spins. The results show exceptional agreement for the final burst of radiation, with some differences attributable to small spins on the black holes in one case.Comment: Revtex 4, 5 pages. Published versio

Spontaneous scalarization with massive fields

We study the effect of a mass term in the spontaneous scalarization of neutron stars, for a wide range of scalar field parameters and neutron star equations of state. Even though massless scalars have been the focus of interest in spontaneous scalarization so far, recent observations of binary systems rule out most of their interesting parameter space. We point out that adding a mass term to the scalar field potential is a natural extension to the model that avoids these observational bounds if the Compton wavelength of the scalar is small compared to the binary separation. Our model is formally similar to the asymmetron scenario recently introduced in application to cosmology, though here we are interested in consequences for neutron stars and thus consider a mass term that does not modify the geometry on cosmological scales. We review the allowed values for the mass and scalarization parameters in the theory given current binary system observations and black hole spin measurements. We show that within the allowed ranges, spontaneous scalarization can have nonperturbative, strong effects that may lead to observable signatures in binary neutron star or black hole-neutron star mergers, or even in isolated neutron stars.This research was supported by NSF Grants No. PHY-1065710 and No. PHY-1305682, NASA Grant No. NNX11AI49G, STFC GR Roller Grant No. ST/L000636/1 (F. M. R.) and the Simons Foundation (F. P.). Computational resources were provided by the Orbital cluster at Princeton University and the COSMOS Shared Memory system at DAMTP, University of Cambridge operated on behalf of the STFC DiRAC HPC Facility. The latter equipment is funded by BIS National E-infrastructure capital Grant No. ST/J005673/1 and STFC Grants No. ST/H008586/1, No. ST/K00333X/1, and No. ST/J001341/1.This is the author accepted manuscript. The final version is available from the American Physical Society via http://dx.doi.org/10.1103/PhysRevD.93.06400

Critical Collapse of the Massless Scalar Field in Axisymmetry

We present results from a numerical study of critical gravitational collapse of axisymmetric distributions of massless scalar field energy. We find threshold behavior that can be described by the spherically symmetric critical solution with axisymmetric perturbations. However, we see indications of a growing, non-spherical mode about the spherically symmetric critical solution. The effect of this instability is that the small asymmetry present in what would otherwise be a spherically symmetric self-similar solution grows. This growth continues until a bifurcation occurs and two distinct regions form on the axis, each resembling the spherically symmetric self-similar solution. The existence of a non-spherical unstable mode is in conflict with previous perturbative results, and we therefore discuss whether such a mode exists in the continuum limit, or whether we are instead seeing a marginally stable mode that is rendered unstable by numerical approximation.Comment: 11 pages, 8 figure

Black Hole Mergers and Unstable Circular Orbits

We describe recent numerical simulations of the merger of a class of equal mass, non-spinning, eccentric binary black hole systems in general relativity. We show that with appropriate fine-tuning of the initial conditions to a region of parameter space we denote the threshold of immediate merger, the binary enters a phase of close interaction in a near-circular orbit, stays there for an amount of time proportional to logarithmic distance from the threshold in parameter space, then either separates or merges to form a single Kerr black hole. To gain a better understanding of this phenomena we study an analogous problem in the evolution of equatorial geodesics about a central Kerr black hole. A similar threshold of capture exists for appropriate classes of initial conditions, and tuning to threshold the geodesics approach one of the unstable circular geodesics of the Kerr spacetime. Remarkably, with a natural mapping of the parameters of the geodesic to that of the equal mass system, the scaling exponent describing the whirl phase of each system turns out to be quite similar. Armed with this lone piece of evidence that an approximate correspondence might exist between near-threshold evolution of geodesics and generic binary mergers, we illustrate how this information can be used to estimate the cross section and energy emitted in the ultra relativistic black hole scattering problem. This could eventually be of use in providing estimates for the related problem of parton collisions at the Large Hadron Collider in extra dimension scenarios where black holes are produced.Comment: 16 pages, 12 figures; updated to coincide with journal versio

Ninja data analysis with a detection pipeline based on the Hilbert-Huang Transform

The Ninja data analysis challenge allowed the study of the sensitivity of data analysis pipelines to binary black hole numerical relativity waveforms in simulated Gaussian noise at the design level of the LIGO observatory and the VIRGO observatory. We analyzed NINJA data with a pipeline based on the Hilbert Huang Transform, utilizing a detection stage and a characterization stage: detection is performed by triggering on excess instantaneous power, characterization is performed by displaying the kernel density enhanced (KD) time-frequency trace of the signal. Using the simulated data based on the two LIGO detectors, we were able to detect 77 signals out of 126 above SNR 5 in coincidence, with 43 missed events characterized by signal to noise ratio SNR less than 10. Characterization of the detected signals revealed the merger part of the waveform in high time and frequency resolution, free from time-frequency uncertainty. We estimated the timelag of the signals between the detectors based on the optimal overlap of the individual KD time-frequency maps, yielding estimates accurate within a fraction of a millisecond for half of the events. A coherent addition of the data sets according to the estimated timelag eventually was used in a characterization of the event.Comment: Accepted for publication in CQG, special issue NRDA proceedings 200